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Inflationary models including vector fields have attracted a great deal of attention over the past 
decade. Such an interest owes to the fact that they might contribute to, or even be fully responsible 
for, the curvature perturbation imprinted in the CMB. However, the necessary breaking of the 
vector field’s conformal invariance during inflation is not without problems. In recent years it 
has been realized that a number of instabilities endangering the consistency of the theory arise 
when the conformal invariance is broken by means of a non-minimal coupling to gravity. In this 
paper we consider a massive vector field non-minimally coupled to gravity through the Gauss- 
Bonnet invariant, and investigate whether the vector can obtain a nearly scale-invariant perturbation 
spectrum while evading the emergence of perturbative instabilities. We find that the strength of the 
coupling must be extremely small if the vector field is to have a chance to contribute to the total 
curvature perturbation. 

PACS numbers: 


I. INTRODUCTION 

Thanks to a wealth of high precision cosmological observations, specially those obtained by the WMAP [H-Q and 
Planck missions M, cosmological inflation is widely recognized as the simplest paradigm to generate the observed 
adiabatic, nearly scale-invariant, Gaussian spectrum of superhorizon fluctuations imprinted in the Cosmic Microwave 
Background (CMB). In particular, single-field models, in which the inflationary expansion is driven by the energy 
density of a scalar field minimally coupled to gravity, are clearly favoured by data. Despite their excellent agreement 
with the available data, indications exist suggesting that single-field models might need to be extended. The most 
notorious among these indications are the so-called CMB anomalies, observed for the first time by the WMAP satellite 
0 and later confirmed by Planck m- 

Although scalar fields have played a dominant role in inflationary cosmology, over the last decade it has been 
realized that vector fields may also have an important function provided their conformal symmetry is broken (see for 
instance [T2| - [T3| and references therein). This breaking, which can be brought about by the introduction of a mass 
term, for example, allows the vector field to obtain a superhorizon spectrum of perturbations during inflation. In turn, 
this opens up the possibility that the vector field becomes a curvaton and contributes to the curvature perturbation 
[l5l - l23l| . for which the vector field must come to dominate (or nearly dominate) the energy density at a later epoch. 
However, the risk when considering the influence of vector fields in the cosmological dynamics is that, since they 
signal a preferred direction in space, they may result in an anisotropic expansion in excess of the current observational 
bounds. To quantify the anisotropy it is usual to parametrize the spectrum of the curvature perturbation as [2^ 

V^{k)=Vl^°ik)[l+g{k){d.kf]^ ( 1 ) 

where V^^°{k) denotes the isotropic part of the power spectrum, d is the unit vector signaling the preferred direction, 

k = k/k is the unit vector along the wavevector k (with k being the modulus of the latter) and g{k) is the so-called 
anisotropy parameter, which quantifies the statistical anisotropy in The latest observations from the WMAP 
and Planck satellites suggest that g can be at most 2% [2^, [^, which represents a very tight constraint on the 
contribution of vector fields to the power spectrum of the CMB. Nevertheless, if the isotropy of the expansion is 
approximately preserved, vector fields could even be responsible for inflation [27H29j. The requirement in this case 
is to have a large number (typically in the hundredths) of randomly oriented vector fields so that, collectively, they 
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result in an isotropic expansion . Remarkably, it is also possible to retain an isotropic inflatio nary ex pan sion while 
dramatically reducing the number of participating vector fields. This is the case of gauge-flation |14l l30l. l3l|. in which 
a non-Abelian gauge field minimally coupled to gravity plays the role of the inflaton. In this proposal, an SU(2) 
gauge field is considered to form a triad of mutually orthogonal vectors, which in turn allows the gauge field to drive 
inflation without giving rise to an anisotropic expansion. 

An interesting manner to break the conformal invariance is by considering a non-minimally coupled vector held, 
thus resulting a modification of gravity [l^ . . Unfortunately, the non-minimal coupling to gravity is known to 

be problematic due to the emergence of instabilities [12, [32 - l33 | . Although the existence of instabilities represents a 
serious drawback for the consistency of the theory, the very nature of the instabilities has been called into question, 
and a number of scenarios have been envisaged to evade them^ [22| . In this paper we examine a cosmological vector 
field non-minimally coupled to gravity through the Gauss-Bonnet invariant. Although couplings between the Gauss- 
Bonnet invariant and scalar fields have been explored in the context of inflationary cosmology [32 - l42| . the coupling 
between the Gauss-Bonnet invariant and a massive vector field has not been explicitly explored in the literature^. 
Arguably, this owes to the very presence of instabilities in relatively simple settings, as in the case of a non-minimal 
coupling to the Ricci scalar, which then invites to exercise caution when considering more complicated non-minimal 
couplings. Nevertheless, the reason for us to invoke such a coupling owes to the peculiar behavior of the Gauss-Bonnet 
invariant. Indeed, a most crucial feature is that it changes sign when passing from inflation to a matter or radiation 
dominated phase. Gonsequently, a mass term for the vector held coming from such a coupling features the same 
change of sign towards the end of inflation. In this regard, the goal of this paper is to investigate whether this change 
of sign allows us to avoid the perturbative instability associated with the change of sign in the vector field’s mass 
squared. Also, we wish to explore whether a vector held with such a coupling can contribute significantly to the total 
energy density after inflation, thus being able to play the role of a curvaton. 

The paper is organized as follows. In Sec.|n]we study the background dynamics of a vector field when coupled to 
the Gauss-Bonnet invariant, and study its perturbation spectrum in Sec. lIIIl In Sec. IIVI we establish the bound on the 
coupling (a) to the vector field necessary so that the latter does not disrupt the inflationary background. In Sec. |V] 
we show that the strength of a must be time-dependent if the vector field is to contribute to the energy density at a 
later epoch. We summarize our conclusions in Sec. ED 


II. VECTOR FIELD COUPLED TO THE GAUSS-BONNET INVARIANT 

We consider a massive vector field coupled to the Gauss-Bonnet invariant and evolving in an inflationary 
background, which we take to be quasi-de Sitter {H ~ 0) and driven by an unspecified matter source. The Lagrangian 
of the system is 

C = Anf - i mlR - + A^A >^, (2) 

where = df_iAi, — d^A^^, R is the Ricci scalar, Q = R^ — ARapR°‘^ + Rap-ysR^^^^ is the Gauss-Bonnet topological 
invariant and a is the strength of its coupling to the vector field A^, whose dimensions are [a] = M~‘^. The effective 
mass squared of the vector field is 

+ aQ , ( 3 ) 

where, in the case of a FRW universe with ds^ = dt^ — a?'{t){dx^ + dy'^ + dz^), the Gauss-Bonnet invariant is 

g = 24:H^{H + H^). ( 4 ) 

Varying the action with respect to the metric we obtain the held equations 

SCeh + + SjCa + = 0, (5) 

where 

SFeh = - ^ {-R^'^ + , ( 6 ) 


^ For other studies regarding the uses and issues of vector fields in cosmology see for example [35l437|| . 

^ Non-minimal couplings of the electromagnetic field to gravity, in particular to the Gauss-Bonnet invariant, have been considered as a 
mechanism to generate large-scale magnetic fields during inflation [4^. [ 4 ^ . 
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5Ca 




( 7 ) 


and 

6Cg = -^agA^^A'^ + i g^^'^fg - 2fRR^'^ + 2V^V'{fR) - 2gf^''D{fR) + 8fRi^R''P - 4VpV'"(/i?‘^'’) 

-AVpV''{fRPP) + AOifRP^^) + 4gP''VpV^{fRP'^) - 2fRPP^^R'^p^^ + AVpS/^{fRP p^'^) , ( 8 ) 

where we define the function / = 5 aA^ = —aA^j2 after choosing a reference frame with the z-axis in the direction 
of the background vector field A. The Friedmann equation that follows can be written as ZH^m^p = pinf + PA^ where 
Pinf is the energy density of the source driving inflation and 


A^ fh?A\ rr4 ^2 / 


2Ha 


41. 


( 9 ) 


is the energy density of the vector field. Note that we have allowed a to be time-dependent. Varying now the action 
with respect to 41^ = (0, A) = (0, 0,0, A.) we obtain the evolution equation for the homogeneous vector field 

A^+HA^+m^A^=Q. ( 10 ) 


In terms of the physical vector field component B = A^ja, the last two equations become 

B + 3HB + (m^ F 2H^ + H)B = 0 . 

and 


1 ■, 1 
Pa-^B +- 


fhr -t- 1 -f 


2B 

llB 


B^ + 2Aaa^B^H‘^ I 1 -1---t 


B a 


HB 2Ha 


( 11 ) 


( 12 ) 


This is to be compared with the energy density pA = +1 rri^B^ of a vector field non-minimally coupled to gravity 

through the Ricci scalar [I^, [27H^ . In our case, the energy density of the vector field crucially depends on the 

scale factor o(t), which can result in a exponential variation during inflation. Owing to this, a salient feature of this 
scenario is that the vector field might be able to retain a significant fraction of its energy density even if it is not in 
slow-roll. In turn, this might be of interest for cosmological purposes. However, as we see below, this dependence 
with the scale factor a{t) prevents the vector field from playing a significant role in a cosmological context. 


III. PERTURBATION SPECTRUM 
A. On the instability of non-minimally coupled vector fields 

An issue of fundamental importance concerning theories of massive vector fields with a non-minimal coupling to 
gravity is that of their instability [13, l32l - [^ . which originates from the longitudinal mode of the vector field. One 
of the known instabilities is perturbative in origin, and arises when the effective mass squared of the vector field 
changes sign from negative to positive In the scenario studied in [13, it was shown in [13 that the 

instability is under control during inflation. However, the instability arises at some later epoch, when the field’s 
effective mass squared crosses zero. In spite of this difficulty, the authors in go on to argue that even if such 
instability exists, it still might be possible to avoid it if the bare mass of the vector field stems from the coupling to 
another field, which then would allow either a curvaton or an inhomogeneous reheating mechanism. To the best of 
our knowledge, the debate on this issue is not yet settled, and hence our attitude towards it will be the same as in 
[^ . thus simply ignoring the instability or assuming that, if present, it can be circumvented by some mechanism. 
Although this attitude conveniently dispenses with the problem, it is also fair to say that such instability arises when 
the longitudinal mode becomes unphysical, and hence it is reasonable to suspect that the associated singularity might 
share the same unphysical nature. 

Apart from the above, yet another problem plagues this kind of vector field models, the so-called ghost instability. 
This originates because, during inflation, the kinetic energy density for the longitudinal modes of the vector field have 
the wrong sign, which might entail the copious production of vector field quanta up to the point of ruining inflation. 
Regarding this instability, the authors in (2^ argue that as long as the negative energy contributed by ghost states 
does not exceed the energy density driving inflation, these are in principle not problematic for the stability of the 
theory. In the following, we implicitly assume that this is indeed the case. 
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B. Perturbation spectrum 


Having clarified our position with respect to the instability of the theory, we investigate the range of a for which the 
vector field obtains a nearly scale-invariant spectrum of superhorizon perturbations. Following |15l |. one can find the 
equations for the temporal (At) and spatial (A) components of the vector field Introducing the physical vector 
field B = A/a, we obtain the equation for B 


SB 3HSB + {2H^ 


- a-^V^) SB + 2H a-^VAt = 0. 


(13) 


In order to quantize the field, we introduce creation/annihilation operators for each polarization 

(Pk 


SBx{t,x) = 


(27r)3 


ex{k)ax{k)bx(t, -b e^(fc)o^(fc)&^(f, k)e 


-ik-x 


(14) 


where k = k/k, k = \k\, A = L,R,\\ labels the left and right transverse and longitudinal polarizations respectively, 
and we use polarization vectors e>(fe) and commutation rules are as in [T^. With all the above, the perturbation 
spectrum for the polarization A is defined by 


Vx{k) = lim 




k/aH^O 2Tr^ 


(15) 


For our purposes, it will be enough to study the transverse components of the vector field. Therefore, we simply 
consider the evolution equation for the mode functions bL,R [l^ 


(16) 


bL,R + -b aO -b 2H -b ^ — 0 , 

whose solution, matching the Bunch-Davies vacuum in the subhorizon limit k/aH oo, is 


bL,R{t, k) = a 


= »-3/2 





where 


= 7 - 
4 


oiQ 




Taking the superhorizon limit k/aH —^ 0 of Eq. (1171) we easily obtain 


Bl.r = 


2-i+2uY{iyy 


k 


3-2!^ 


The spectral index of the field perturbation spectrum is defined by 

nL,R - 1 = 3 - 2i/, 


(17) 


(18) 


(19) 


( 20 ) 


with scale-invariance corresponding to ul^r = 1. From Eqs. (ITOl) and (fTOl) . and according to the findings in [T^, it 
is clear that only when ss —2H^ can the vector field be substantially produced during inflation while attaining a 
nearly flat perturbation spectrum. Consequently, the coupling a must satisfy -b aQ ~ —2H^. As is well-known 
when a = 0, a flat perturbation spectrum can be obtained by coupling the vector field to the Ricci scalar [Il,[2i,[23- 

S However, as previously discussed, this scenario is known to have a number of problems concerning its stability 
l32l - [^ . A natural question to investigate then, is whether the coupling to the Gauss-Bonnet invariant can alleviate 
such problems to some extent. 


IV. INCONSISTENCY OF A SCALE-INVARIANT SPECTRUM 

From Eq. m it follows that the condition m? ~ —2H^ to achieve a flat perturbation spectrum implies a slow-roll 
motion for the homogeneous field B. In this case, the contribution of the Gauss-Bonnet term to the energy density of 
the vector field, given by the last term in Eq. (HU, grows exponentially during inflation, provided a is a constant. As 
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a result, unless a is extremely small, the natural outcome is that the rapidly growing energy density pA will disrupt 
the phase of slow-roll inflation. Therefore, in order to keep the inflationary background oblivious of the vector field 
it is necessary to control the growth of pA- In view of Eq. (HID, the simplest alternative is to enforce the condition 
aa^B^ ~ const. In such case, if a remains constant during inflation the homogeneous field B cannot be in slow-roll, 
but it must perform a fast-roll motion. To be precise, writing B oc a~^ the condition a^B^ = const, becomes ,5 ~ 1. 
Let us now define S by a^B^ oc (namely S = 2 — 2/3), so that the condition /3 ~ 1 can be rewritten as |<3| < 1, 
and also let Ntot be the total length of inflation. Using these, the field’s energy density pA grows by a factor of order 
exp(3Vtot^) by the end of inflation, whereas the first two terms in Eq. (I12p . while proportional to B^, decay by a factor 
of order exp(—2A^tot), and hence we can neglect them. As a result, we can approximate the field’s energy density by 

PA ^ l2{a5)B^H‘^ exp(iV,5), (21) 

where Bi is the value of B at the beginning of inflation and N is the number of e-foldings. From this expression we 
see that in order to have > 0 we must enforce the condition 0(5 > 0. If IVtot is not substantially larger than the 
number of e-foldings necessary to comply with observations (typically around 60) and 5 is relatively small, then pA 
may remain subdominant during the entire inflationary phase. Demanding that pA < Pinf — where iL, is 

the Hubble scale during inflation, holds until the end of inflation we obtain 

(“^')^exp(3V,o,5)<l(|i) . ( 22 ) 

To relate the above constraint with the perturbation spectrum, we use the growing mode of the solution to Eq. m- 
With B (X a ^, this is determined by 


-P = 


3 

2 



ifi^ + aG 
IP 



(23) 


and since (5 = 2 — 2/3, we can use Eqs. (ITSl) and (l20)) to find S — 2 — — 1). For a flat perturbation spectrum 

((5 ~ 2 ), and for fields able to be produced during inflation (for which a natural assumption is that B'^ ~ the 
constraint in Eq. (1221) implies that a must be extremely small. Such a constraint is significantly alleviated when iVtot 
is set to its minimum at Wot — W and when \Bi\ <C H^, at the onset of inflation. Nevertheless, the level of tuning of 
a is still quite large. However, the most effective manner to reduce the tuning is to consider |(5| <C 1, which implies a 
strongly scale-dependent spectrum for B. Below we study this case in detail. 

Firstly we examine the case vP > 0. From Eq. (l23)) and the definition of 8 we have 


vP -I- aG 





(24) 


where the last approximation follows from |(5| <C 1. Using now the condition > 0 (necessary to have pA > 0) in the 
above we obtain |a|t3 < |to^|, and hence rP > 0. This case is desirable if one is to avoid the perturbative instability, 
for rP does not change sign at any time during the entire evolution. In principle, for tP > 0 and a < 0 the constraint 
in Eq. (l2^ leaves plenty space available (which coincides with the shaded region in Fig. [T]). However, from Eq. (l24ll 
we see that > 0 corresponds to (5 < 0, hence pA decreases during inflation (see Eq. (HU). Also, in view of Eqs. (fT51) 
and ((^ we obtain > 3. Using for example Wot = 100 and W = 10“®top and restricting ourselves to the range 
10“® < aHl < 1 we find 3 < np. R ^ 4. As a direct result of the evolution implied by the strong scale dependence and 
the decrease of pA, the field’s energy density by the end of inflation is too much suppressed (as anticipated), which 
then makes impossible that the vector held contributes to the total energy density at any later epoch. 

Our results for fp < 0 and q( > 0 are illustrated in Fig. [TJ Since fP < 0 corresponds to (5 > 0, the held’s energy 
density pA grows during inflation. The shape of the shaded region, corresponding to the range of parameters allowed 
by the constraint in Eq. (1^^ . is very sensitive to Wot due to the exponential factor. For the case shown in Fig. [TJ 
taking Wot = 100 and W = 10“®top (corresponding to ^ 10^'^GeV) we hnd 2.5 < ul^r < 3. The right-hand 
panel evidences that a must be extremely small if the vector held is to obtain a hat perturbation spectrum, which 
corresponds to fP = —2H^. The region with heavy shading is where the vector density parameter flA becomes of 
order 1 at some later epoch after inhation, thus giving a sizable contribution to the total energy density (see Sec. IVl). 


V. THE NEED FOR A TIME-DEPENDENT a 

As seen in the previous section, during inhation pA is dominated by the Gauss-Bonnet contribution. Since pA oc 
its scaling after inhation is very fast. Indeed, if a remains constant then pA oc a~^B^ if the Universe becomes radiation 
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FIG. 1: The figure depicts the range of masses < 0 and couplings a > 0 for which pA remains subdominant during the 
entire inflationary phase (shaded area). To build the plot we used A^tot = 100 and IT, = 10~®mp. In the left-hand panel, the 
spectral index in the shaded region is 2.5 < nL,R < 3. The region with heavy shading represents the case when the vector field 
manages to give a sizable contribution to the total energy density at some later epoch. The right-hand panel evidences the 
extreme smallness of the coupling a necessary for the vector field to obtain a scale-invariant perturbation spectrum. 


dominated after inflation, which demands a rapid growth of the homogeneous field B if the latter is to give a sizable 
contribution to the total energy density at a later epoch. Such a stage of rapid growth is difficult to motivate, although 
one could envisage the case of a vector field with a symmetry breaking potential, as is the case of bumblebee models 
[dsl - f^ . In these, at least in principle, one might arrange that the field B grows towards its vacuum expectation value 
during the post-inflationary epoch. We do not pursue this alternative here. Instead, we investigate the case when the 
coupling a is time-dependent. A possible motivation for this is to conjecture that the magnitude of a is determined 
by the effective vacuum energy density, thus remaining constant during inflation and scaling as a oc H~^ oc if 
the Universe becomes dominated by the matter-like oscillations of the inflaton. Below we focus on this case, thus 
assuming that 


aH^ = const., (25) 

but without restricting our analysis to a particular physical motivation. This is to say that, in the interest of generality, 
we maintain the equation of state parameter of the fluid w that dominates after inflation. 

After inflation, when the energy density is dominated by a fluid with equation of state parameter w, the energy 
density of the vector held reads 


PA = 2Aaa^B^H^ (^1 + 7- + , (26) 

where we have written B oc a'''™ as the formal solution to Eq. (HU. After inflation, the contribution of the Gauss- 
Bonnet term to the field’s mass squared scales as aQ oc . Moreover, since due to the condition > 0 the mass 
squared term aQ cannot dominate m^, given the rapid scaling of aQ after inflation we can approximate m? ~ fri^ and 
Eq. ((TT]) by 

B + ?,HB+(^'^ + ^H^{l-iw)^B = Q. (27) 

When ^ the field starts performing fast oscillations around the minimum of its potential, hence = —3/2 
and Pa/oc. a^B^ oc. a~^. Asa result, the vector field never gives a significant contribution to the total energy density. 
However, before the vector field engages into the oscillatory regime, i.e. when we find = —(I — 3w)/2. 

Therefore, the vector field has a chance to contribute to the total energy density provided QIa = Pa/ P oc 
grows with time. This happens for w > —1/3, which encompasses the case when the Universe becomes dominated by 
the matter-like oscillations of the inflaton. Taking into account that oc and demanding that the vector 
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field gives a sizable contribution to the total energy when H = Hgq, i.e. imposing ~ 1, the density parameter 

VIa at the end of inflation is constrained by 


nA,e 4(ai7^) ^ exp {Ntot^) 


m) 


4(1+7x(;) 

3{l+in) 


(28) 


where we used Eq. m to approximate 12^,6- Of course, the vector field manages to contribute to the total energy 
only if riA keeps growing for a sufficiently long-lasting period, which is obviously guaranteed if 


< H, 


eq 


^2 (2(l-|-7m)) 


(29) 


This condition, however, is weaker than the constraint in Eq. (1221) and hence it does not result in a reduction of the 
available parameter space. In the case shown in Fig. (TJ the region with heavy shading corresponds to a vector field 
giving a sizable contribution to the total energy density when iJgq ^ GeV, well before the time of Nucleosynthesis, 
which is the latest time when a matter field with a sizable fraction of the total energy can decay into radiation. 


VI. CONCLUSIONS 

In this paper we have examined the evolution of a cosmological vector field coupled to the Gauss-Bonnet invariant. 
We find that in the presence of such coupling, the energy density of the vector field grows exponentially during inflation 
unless B decreases approximately as a~^. In that case, however, the vector field obtains a strongly scale-dependent 
perturbation spectrum. A flat perturbation spectrum is obtained only at the expense of a severe tuning of the coupling 
a and a limited, total number of e-folds during inflation (right-hand panel in Fig.[T]). We also find that if the strength 
of the coupling a remains constant after inflation, the field’s energy density decreases too quickly to contribute to the 
total energy density at any later epoch. In order to allow the vector field to play a role at a later time, we assume 
a time-dependent coupling scaling as the inverse of H^, which at the same time implies a nearly constant a during 
inflation. We show that if a is sufficiently small, the vector field may come to dominate the Universe at a later time 
(see Fig.[T|). 

As a final comment, we note that the flatness of the perturbation spectrum is a necessary requirement when 
the vector field is responsible for the most part of the curvature perturbation. Nevertheless, as shown in , a non- 
minimally coupled vector field cannot do so, for it gives rise to a spectrum too anisotropic, with g ~ — 1. Gonsequently, 
although the requirement of a nearly scale-invariant perturbation spectrum can be relaxed, the scale dependence found 
for the space depicted in Fig. [1] is too large to be of interest. 
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